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Surface-confined 2D polymerization of a brominated copper-tetra- 
phenylporphyrin on Au(lll) 

Lars Smykalla,*" Pavel Shukrynau," Marcus Korb,* Heinrich Lang,^ and Michael Hietschold" 


A coupling-limited approach for the Ullmann reaction-like on-surface synthesis of a two-dimensional covalent organic network 
starting from a halogenated metallo-porphyrin is demonstrated. Copper-octabromo-tetraphenylporphyrin molecules can diffuse 
and self-assemble when adsorbed on the inert Au(lll) surface. Splitting-off of bromine atoms bonded at the macrocyclic 
core of the porphyrin starts at room temperature after the deposition and is monitored by X-ray photoelectron spectroscopy 
for different annealing steps. Direct coupling between the reactive carbon sites of the molecules is, however, hindered by the 
molecular shape. This leads initially to an ordered non-covalently interconnected supramolecular structure. Further heating to 
300 °C and an additional hydrogen dissociation step is required to link the molecular macrocycles via a phenyl group and form 
large ordered polymeric networks. This approach leads to a close-packed covalently bonded network of overall good quality. 
The structures are characterized using scanning tunneling microscopy. Different kinds of lattice defects and, furthermore, the 
impact of polymerization on the HOMO-LUMO gap are discussed. Density functional theory calculations corroborate the 
interpretations and give further insight into the adsorption of the debrominated molecule on the surface and the geometry and 
coupling reaction of the polymeric structure. 


1 Introduction 

The formation of molecular nano-structures on surfaces 
through self-assembly of building blocks which are held to¬ 
gether only by weak intermolecular interactions has the in¬ 
herent problems that these arrangements are fragile and in¬ 
termolecular charge transport is limited. A solution to this is 
a controlled interconnection of the molecules on the surface 
by robust and irreversible covalent bonds. Covalent linking 
of small organic molecules directly on the surface allows also 
for the engineering of manifold extended 2D materials which 
are otherwise not accessible by classical solution-based poly¬ 
merization. Especially, the prospect to fabricate ;r-conjugated 
graphene-like nano-sheets of controlled dimensions is of great 
interest. Therefore, the on-surface synthesis of ID and 2D 
polymers has attracted considerable attention over the recent 
years. ^ ® 

One major approach for surface-confined polymerization is 
the halogen-based covalent self-assembly based on the Ull¬ 
mann coupling reaction. ^ Halogen atom substituents are split- 
off a monomer and, thus, unsaturated carbon atoms at pre¬ 
defined positions are created. Thereby, the influence of a 
catalyst (originally Cu) or the adsorption on a surface sig¬ 
nificantly lowers the energy barrier for dissociation. Subse¬ 
quently, the resulting molecular radicals diffuse to fitting po¬ 
sitions close to another to enable C-C coupling into covalent 
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dimers, chains or networks. The topography and also elec¬ 
tronic properties of the covalent organic framework are de¬ 
termined by the monomer used and its halogen substitution 
pattern. Also, the properties of the substrate surface, like the 
reactivity and the energy barriers for diffusion and rotation, 
play an important role in this kind of synthesis. 

The polymerization is often observed to not occur directly 
after the dehalogenation, rather an additional activation step 
is required, which is achieved typically photochemically or 
by heating. Also, it was demonstrated that the reaction steps, 
which are halogen dissociation, movement of the molecules 
and C-C coupling, could even be individually and locally per¬ 
formed by voltage pulses and manipulation with the tip of a 
scanning tunneling microscope. 

It was reported that dehalogenation started already at re¬ 
duced temperature and is nearly finished at room temper¬ 
ature especially on the reactive Cu surfaces. After¬ 
wards, the radicals can be stabilized in organometallic struc¬ 
tures by metal atoms, e.g. Cu surface adatoms. 

For actual polymerization, additional heating often to high 
temperatures is then necessary to cleave the carbon-metal 
bonds and enable diffusion that subsequently leads to irre¬ 
versible C-C coupling of the dehalogenated radicals. This 

diffusion-limited polymerization route, however, leads to the 
typical main problems of 2D polymers, which are small do¬ 
main sizes a high density of defects or undesired crum¬ 
pled or dendritic networks without long-range order. 
Split-off halogen atoms which remain on the surface stabi¬ 
lized by chemisorption and the presence of metal adatoms 
can also hinder the formation of long-range ordered covalent 


1 





networks Considering that a high mobility of molecules 
is a prerequisite for the formation of well-ordered structures, 
Cu(lll) seems also not to be the ideal surface for halogen- 
based covalent self-assembly due to its relatively high dif¬ 
fusion barriers. The optimal conditions for large ordered 
2D networks are a surface with lower diffusion barriers like 
Au(lll) as well as a relatively large recombination barrier 
of the monomer radicals, which results in a coupling-limited 
polymerization process. This means, monomers are able 
to diffuse and rearrange to a greater extent before they finally 
irreversibly combine into polymeric networks. Notably, prob¬ 
ably the first example of halogen-based, two-dimensionally 
covalent self-assembly on a surface was demonstrated on 
Au(lll) by Grill and coworkers.^^ However, higher temper¬ 
atures to effectively initiate the reaction steps are necessary 
on Au compared to Cu surfaces, which could already lead to 
further decomposition of large monomers. 

In general, scanning tunneling microscopy (STM) has been 
the main tool for exploring surface-confined polymerization. 
Although, despite the formation and structure of various 2D 
polymers were reported, only little attention was given to 
their electronic structure.^ Compared to covalently bonded 
networks of metal-free molecules, even more interesting trans¬ 
port properties are expected, when incorporated metal ions 
can be coupled. For applications in molecular spintronic de¬ 
vices, the spin bearing metal ions should be decoupled from 
the surface, which mostly excludes the metal-organic coor¬ 
dination assemblies linked by metal adatoms. A good 
candidate could be networks of metallo-tetraphenylporphyrin 
molecules which additionally should have a larger mechan¬ 
ical robustness and stability on the surface compared to non- 
covalently bonded structures of “magnetic” molecules. Nev¬ 
ertheless, the largest challenge currently still lies in improv¬ 
ing the domain size and quality of the polymeric organic net¬ 
works. Previously, Grill and coworkers investigated a metal- 
free tetraphenylporphyrin with bromine at the para position 
of the phenyl groups which couples the molecules at these 
C atoms of the phenyl rings after the cleavage of the C-Br 
bonds and, consequently, forms networks with a low pack¬ 
ing density. In this paper, we report the steps for polymer¬ 
ization of copper-2,3,7,8,12,13,17,18-octabromo-5,10,15,20- 
tetraphenylporphyrin (CuTPPBrg) monitored by a combina¬ 
tion of STM and X-ray photoelectron spectroscopy (XPS) and 
supported by density functional theory calculations. The dif¬ 
ferent bromine substitution results in a hindered direct C-C 
coupling and, thus, a new coupling scheme and close-packed 
covalently bonded network, which will be discussed in the fol¬ 
lowing. 



Fig. 1 Chemical structure of copper-2,3,7,8,12,13,17,18- 
octabromo-5,10,15,20-tetraphenylporphyrin 
(CuTPPBrg) 

2 Experimental details 

A clean surface of the single crystal of Au(l 11) was prepared 
by multiple cycles of Ar+-sputtering at an energy of 500 eV 
and annealing to 400 °C for 1 h. CuTPPBrg was synthesized 
by a similar procedure as reported in Ref. 29. The chemical 
structure is shown in Fig. 1. The molecules were deposited 
on the substrate by organic molecular beam epitaxy in ultra 
high vacuum (UHV). Before, the evaporant was first purified 
by heating to a temperature slightly below the sublimation 
temperature in UHV until the pressure dropped again. The 
molecules were then deposited at around 350°C on Au(lll). 
The substrate was kept near room temperature during deposi¬ 
tion. Scanning tunneling microscopy experiments were per¬ 
formed with a variable temperature STM from Omicron in 
UHV. The base pressure in the UHV chamber was in the 
range of 10“^^mbar during measurements. Electrochemically 
etched tungsten tips for the STM were cleaned by annealing in 
UHV. All measurements were done at room temperature and 
STM images were recorded in constant current mode with a 
tunneling current of lOOpA if not specified otherwise. Pos¬ 
itive bias corresponds to electron tunneling from the tip into 
empty states of the sample. STM images were processed with 
the WSxM software. 

X-ray photoelectron spectroscopy experiments were per¬ 
formed with synchrotron radiation at the Material Science 
beam line in Elettra (Trieste, Italy). A Phoibos photoelectron 
spectrometer was used. The excitation photon energy was for 
Br 3d 180eV, for C Is 400eV, and for N Is 520eV. The bind¬ 
ing energies were corrected relative to the Au 4f7/2 peak of the 
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substrate located at 84 eV. The shapes of the core-level signals 
were fitted using Voigt peaks and a Shirley background. The 
samples were heated up to the temperature indicated in the 
figures, which was measured directly at the backside of the 
crystal, held at this temperature for 10 min and then cooled 
down to room temperature prior to the XPS measurement. 

Density Functional Theory (DFT) calculations were carried 
out using the grid-based projector augmented wave method. 
The RPBE exchange-correlation functional with a pair-wise 
correction for the dispersion interaction [vdW(TS)] and the 
vdW-parameters for Au by Ruiz et al were used to include 
screening effects of the metal surface. A grid spacing of 
h = 0.2 A, and a dzp basis set of atomic orbitals for the wave 
functions were used. For the self-interaction of the d states of 
the Cu atoms a Hubbard U correction with a value off/ = 5 eV 
was applied. Simulation cells with periodic boundary condi¬ 
tions in X and y and zero boundary in z direction with a total 
vacuum of both sides of 1.8nm were used. For the calcula¬ 
tion of debrominated CuTPP on Au(l 11), a slab of 3 layers of 
Au was used, whereby the two lowest layers were fixed during 
optimization. The reconstruction of the Au(lll) surface was 
not considered within the calculations. The structures were re¬ 
laxed until all atomic forces converged. Only the F-point was 
used for the Brillouin zone sampling due to the relative large 
size of the cell. 

3 Results and discussion 

3.1 XPS investigations 

The evolution of the XPS core-level spectra with increasing 
temperature is shown in Fig 2. For Br 3d, the shape has to 
be fitted with three components if an identical width and spin- 
orbit splitting (1.04eV) is used for the doublet peaks. The 
doublets with the larger Br 3d5/2 peaks located at 67.7 eV and 

68.1 eV (black and pink curves) are attributed to Br 2 ^^ and 

Br chemisorbed on the surface. 16,18,19,35 doublet at the 
higher binding energy with the Br 3d5 /2 peak at 69.7 eV (green 
curve) was reported to correspond to bromine bonded to car¬ 
bon 16,19,34, I ^ gj. intact CuTPPBrg molecule. 

After deposition at room temperature, a high amount of disso¬ 
ciated bromine (doublets at lower binding energy) was found. 
It is to be expected that the Br atoms in CuTPPBrg which are 
in direct contact to the surface have a lower energy barrier for 
dissociation than the Br atoms which are bended away from 
the surface due to the saddle deformation of the molecule. It 
is also likely that some molecules already partially debromi¬ 
nated during the sublimation in the crucible or when arriving 
with high thermal energy on the sample surface. The dis¬ 
sociated bromine atoms, which are chemisorbed on the Au 
surface, could then also recombine with the radical molecules 
or form Br 2 . 
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Fig. 2 Evolution of the XPS spectra of Br 3d (a), C Is (b), and N Is 
(c) for a molecular monolayer of CuTPPBrg on Au(lll) compared 
between annealing steps to the indicated maximal temperature. In 
the Br 3d spectra, the background is subtracted. 


After heating the sample to 130°C for 10 min, the C-Br 
peaks were reduced while the other doublets slightly increased 
in intensity, which means that further molecules were de¬ 
brominated. If the sample is annealed to a maximum of 
200 °C, the debromination of all molecules is finished and 
Br atoms or Br 2 remain adsorbed on the surface but also al¬ 
ready start to desorb from Au(lll). Finally, after heating to 
300 °C, no Br 3d signal was found anymore, i.e. bromine com¬ 
pletely disappeared from the Au(lll) surface. Furthermore, 
partial desorption (coverage reduced by 25%) after annealing 
to 300 °C is indicated by the increase of intensity of the XPS 
peaks of the Au substrate (Au 4f) which are attenuated by the 
molecular layer. The temperature necessary for debromination 
on the relatively inert Au(lll) is higher than those reported 
for other brominated molecules on Cu or Ag surfaces, where 
cleavage of the C-Br bonds was found to be complete at or 
even below room temperature, or at slightly elevated tempera¬ 
tures, respectively. In studies on Au(lll), debromination 
and C-C coupling for a tetraphenylporphyrin with a different 
bromination scheme than CuTPPBrg was, however, found in 
STM studies to occur only after heating to 300 

The heating induced debromination can also be observed in 
the C Is core level spectra [Fig. 2(b)]. A peak component at 
« 284.5 eV (arrows) shrinks with increasing temperature. We 
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attribute this feature to carbon bonded to bromine. However, 
this species overlaps with the peak for C-N at almost the same 
energy, which makes a non-ambiguous deconvolution of com¬ 
ponent intensities difficult. The shift of the C Is (and also N 
Is) signal to higher binding energy by 0.2eV after heating to 
300 °C is related to a shift of the Fermi level due to the des¬ 
orption of bromine from the surface and the polymerization of 
the molecules, which would cancel a possibly strong interac¬ 
tion of unsaturated carbon atoms with the surface. Annealing 
to 400 °C leads to no further changes in the XPS signals. 

Small peaks at higher binding energies compared to the 
main core-level peaks are often shake-up satellites due to an 
additional n ^ 7t^ excitation. For N Is [Fig. 2(c)] with only 
a single component (N-Cu), the satellite can be clearly iden¬ 
tified. Because the N orbitals contribute to the highest occu¬ 
pied (HOMO) as well as lowest unoccupied (LUMO) molec¬ 
ular orbital, the excitation energy can be correlated with the 
energy gap between the HOMO and LUMO.^^ However, the 
distance of the shake-up satellite to the corresponding core¬ 
level peak could be smaller than the actual HOMO-LUMO 
gap due to reorganization of the charge within the molecule 
upon photoexcitation and, thus, changed final state screen¬ 
ing of the core hole.^^ The separation from the main N Is 
peak at 398 eV is 1.2eV for the completely debrominated 
molecules. Heating to 300 °C leads, besides the -1-0.2 eV shift 
of all peaks related to the molecule, to a larger shift of the 
satellite peak indicating a higher HOMO ^ LUMO excita¬ 
tion energy of 1.6eV. This means that by polymerization the 
band gap widened by ^0.4 eV. Gas phase calculations us¬ 
ing DFT [RPBE-FU(5eV)-FvdW(TS)] give a HOMO-LUMO 
gap of 1.45 eV for CuTPPBrg, 0.86eV for the debrominated 
molecule CuTPP, L04eV for an unsaturated CuTPP-polymer 
and 1.58eV for the H-saturated CuTPP-polymer. It should be 
noted that DFT generally underestimates energy gaps. On the 
other hand, the HOMO-LUMO gap is significantly reduced 
in the proximity of a surface compared to a thick layer or 
the gas phase as calculated here, which leads to the seem¬ 
ingly good agreement for CuTPP and the H-saturated CuTPP- 
polymer with the experimental values. After debromination 
the HOMO-LUMO gap is strongly reduced due to C atoms 
with unsaturated bonds. After C-C coupling of the CuTPP 
monomers the gap is increased by circa 0.2 eV if the C atoms 
in the polymer would still be unsaturated. A value slightly 
higher than the energy gap of a single CuTPPBrg is found if 
the pyrrolic C atoms of the final polymer are fully saturated 
with the H atoms previously bonded at the phenyl groups. In 
recent studies, a trend of band-gap contraction with increasing 
sheet size of various 2D polymers was reported in theoreti¬ 
cal gas phase calculations and scanning tunneling spec¬ 
troscopy"^^ due to an extended ;r-conjugation. We propose 
that the observed opening of the HOMO-LUMO gap after C- 
C coupling is also due to the reduced hybridization of orbitals 



Fig. 3 (a) STM image after deposition of^l ML of CuTPPBrg and 
annealing to max. 200 °C. Two periodically ordered rotational 
domains and disordered areas in between are shown. The 
herring-bone reconstruction of the underlying Au(lll) is visible 
through the large ordered areas. (U = —0.7 V) (b) Magnification of 
the adlayer structure of CuTPPBrg on Au(lll) with elongated 
appearance of molecules at t/ = — 1V. The cyan arrow indicates the 
direction of an Au lattice vector, (c) Model of the molecular 
arrangement; Protruding C atoms of the saddle deformation are 
colored yellow. 

of the unsaturated carbon atoms with the Au surface. Fur¬ 
thermore, unlike completely ;r-conjugated networks such as 
graphene, the electron delocalization in the described CuTPP 
polymer is likely limited because of a large twist angle be¬ 
tween the p-orbitals in the phenyl groups and in the porphyrin 
macrocycle, which separates the respective ;r-systems. 

3.2 STM investigations of the molecular adlayer 

In the following, the observations in STM of the changes 
in the molecular arrangement after annealing will be dis¬ 
cussed. After deposition of a coverage of around one mono- 
layer of CuTPPBrg on Au(lll), many disordered areas but 
also small, periodically structured, close packed arrangements 
were found. Annealing of the sample drastically enlarged such 
ordered areas due to the increased mobility of the molecules, 
which enables vast diffusion, rotation and, therewith, rear¬ 
rangement and structure formation. After heating to 200 °C, 
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Fig. 4 (a) New structure of the alternating molecular rows formed after annealing to 350 °C {U = —1.5 V). (b) Overlay of molecular model on 
the structure {U = +1.5 V, 6.5 nm x 6.5 nm). (c) Image of the occupied molecular orbitals at the left side of the ordered structure shown in (a). 
The green circle marks a lattice defect, (d) Defect along molecular rows (marked blue) from the incomplete C-C coupling between all 
molecules, (e) Model of the arrangement with covalent bonds between the debrominated molecules (red circle). Incomplete polymerization is 
marked by the green circle. Dissociated hydrogen atoms from the phenyl groups could diffuse and bond to the still unsaturated carbon atoms 
(arrows). Upwards directed C atoms in the pyrrole units of the porphyrin saddle-shape are colored in bright and dark yellow to indicate 
slightly different angles as explain in the text. 


most of the surface area was covered with the self-assembled 
molecular pattern, however, with still disordered regions be¬ 
tween ordered domains. Fig. 3(a) shows two domains of 
the molecular structure which are rotated by 120° relative to 
each other due to the symmetry of the underlying Au(lll) 
substrate. The "herring-bone" pattern of the 22 x \/3 recon¬ 
struction of Au(lll) was visible through the molecular layer. 
This enables the direct measurement of the rotation of the 
molecular lattice relative to a substrate lattice vector, which 
is found to be e(A) = (45 ±5)°. At a bias ofU = —1V (oc¬ 
cupied states), the molecules appeared in an elongated shape 
[Fig. 3(b)], which reveals that the orientation of the molecules 
is rotated alternately by 90° in a checker-board-like man¬ 
ner. Such an appearance is characteristic for various metaho- 
tetraphenylporphyrin (TPP) molecules due to electron density 
at the site of the central metal ion and the well-known saddle- 
like deformation of the porphyrin macrocycle when adsorbed 
on surfaces.The distances between the molecules were 
measured to he A = B = (I.3± O.I) nm, the angle of the unit 
cell a = (88 ± 2)° and the angle between a lattice vector and 
the orientation of the molecular saddle shape 5 = (10±5)°. 
The model of the supramolecular arrangement derived from 
these values is shown in Fig. 3(c). The arrangement is dom¬ 
inated by the attractive interaction of the edge-to-face conhg- 
uration of adjacent phenyl groupswhich is typical for the 
self-assembly of various TPP on noble metal surfaces. 
However, in these reported structures, the azimuthal orienta¬ 
tion is usually the same for ah TPP molecules which is not the 


case for the debrominated CuTPP molecules in the arrange¬ 
ment described here. 

The XPS results indicate that after heating to 200 °C, ah 
CuTPPBrg molecules should have been debrominated and 
thus be surface-stabilized radicals. Despite that, the molecules 
didn’t polymerize but formed a non-covalent structure (Fig. 3) 
instead. This behavior can be understood by taking into ac¬ 
count that, due to the molecular geometry, a C-C coupling be¬ 
tween unsaturated pyrrolic carbon of neighboring molecules 
is strongly hindered by the phenyl groups. In this sense, af¬ 
ter cleaving of the C-Br bonds, the reactive sites are pro¬ 
tected against immediate polymerization because of steric rea¬ 
sons. The monomers can, thus, still freely diffuse and rear¬ 
range, which is highly desired for the formation of long range- 
ordered, close-packed networks with a small amount of de¬ 
fects. 

After the sample was annealed to around 350 °C for Ih, 
STM imaging revealed that the molecular lattice completely 
transformed into a new structure, which is shown in Fig 4. 
This new molecular network consists of apparently alternating 
rows of molecules with slightly different appearance, thus, ad¬ 
jacent molecules are not equivalent anymore. The distance be¬ 
tween neighboring molecules is similar to the previously dis¬ 
cussed structure in Fig. 3. However, ah molecules are rotated 
relative to each other when compared to the arrangement ob¬ 
served for annealing to max. 200 °C - the elongated shapes at 
U = —1.5 V in one type of row are now oriented along the di¬ 
rection of the row [Fig. 4(a)]. The model of the newly formed 
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molecular structure according to the measured distances and 
angles is shown in Fig. 4(e) and overlayed on the STM im¬ 
age in Fig. 4(b). A 90° rotation of adjacent molecules avoids 
cross-like overlapping of phenyl groups and results in favor¬ 
able parallel-displaced n-n stacking. The phenyl groups 
of adjacent molecules point to and are very close to the C 
atoms of the macrocycle where the bond to Br was cleaved. 
Thus, intermolecular covalent C-C bonds between an unsatu¬ 
rated pyrrolic C atom and a phenyl group were formed and the 
molecules were linked together. To make this possible, one 
C-H bond per phenyl group has to be cleaved. It is very 
likely that this dissociated hydrogen atom has been transfered 
inside the molecule along the phenyl group to one of the still 
unsaturated C sites as illustrated in Fig. 4(e), which will later 
be discussed further. The different appearance at the same bias 
voltage of every second molecule (“rows”) can be explained 
by the alternating orientation of the deformed macrocycle and 
respective bonding of phenyl to pyrrole units bended toward 
or away from the Au(l 11) surface. For every second molecule 
the pyrrole rings close to the substrate have to be bended up or 
the protruding pyrrole units slightly downwards toward the av¬ 
erage molecular plane to minimize the stress on the new C-C 
bond with a phenyl ring of another molecule. Thus, the re¬ 
sulting saddle-shape deformations are non-equivalent for ad¬ 
jacent molecules, which was observed in STM [Fig. 4(a)]. A 
fusion of the pyrrolic C with the phenyl groups of the same 
molecule can be excluded in this case because this would 
result in the complete disappearance of a saddle-shape. 

Each polymeric network covers a relatively large area. 
Nonetheless, coupling defects, which are typical for cova¬ 
lent networks, are still a big problem and will, thus, be dis¬ 
cussed in the following. Four types of structural defects can 
be found. First, there are missing molecules creating holes in 
the network, which also occurs in the non-covalently bonded 
arrangement. On the lower left end of the ordered molecu¬ 
lar structure in Fig. 4(a), two rows with molecules showing 
the elongated shape are next to each other and also not dis¬ 
placed along the direction of the row as usual. The STM im¬ 
age of the unoccupied orbitals [Fig. 4(c)] shows a high den¬ 
sity of states between each pair of neighboring molecules, 
which is an indication for a C-C coupling between the unsatu¬ 
rated carbon atoms of two pyrrole groups. In the same way, 
two molecules inside the same row [green circle in Fig. 4(c)] 
can be found closer together and link directly both porphyrin 
macrocycles. However, as discussed previously, for this a 
strong, energetically unfavorable deformation by bending the 
phenyl groups away would be necessary. The third type of de¬ 
fect, which occurs the most often, is shown in Fig. 4(d) and 
can also be seen magnified in Fig. 4(b): The molecules in ev¬ 
ery second row can have different shapes which are related 
to bonding with all phenyl groups to neighboring molecules 
(“close-packed” appearance) or C-C coupling with only two 


phenyl groups (slightly displaced and more elongated molec¬ 
ular shape). Defects should be problematic for applications 
in devices using lateral electronic transport, but at least the 
last type should be possible to cure by further annealing 
to promote complete C-C coupling between all neighboring 
molecules. 

3.3 DFT calculations 

In this section, we present DFT calculations [RPBE ^-U{5 eV) 

vdW(TS)] of the debrominated CuTPPBrg molecules and 
of the polymeric structure. In Fig. 5, the optimized geome¬ 
try is shown of a debrominated CuTPPBrg molecule (CuTPP) 
adsorbed on the favorable position with the Cu atom of the 
molecule 2.8 A above a bridge position of Au(lll). After 
the relaxation, the molecule has average angles of the phenyl 
planes of 0ph ~ 18°, and of the pyrrole planes of 0pyr ~ 19° 
and ~ 37° relative to the surface plane. Two pyrrole 
rings are, thus, strongly bent toward the surface to create C- 
Au bonds which stabilize the radical molecule, while the other 
pair of pyrrole rings remain with unsaturated carbon atoms. 
The strong saddle-shape deformation is responsible for the 
elongated appearance in Fig. 3. The deformation energy to 
lift up the pyrrole rings which are bonded to the Au surface 
was calculated to be 2.6eV. The necessity of an additional 
heating step for the polymerization could be due to the break¬ 
ing of C-Au bonds and reversal of the deformation of the pyr¬ 
role rings and, thus, reactivation of these unsaturated carbon 
atoms. Furthermore, rotational and diffusion energy barriers 
have to be overcome for the reorganization of the molecular 
arrangement. Also, a high annealing temperature could be 
required due to bromine atoms which could be located be¬ 
tween the molecules after debromination (although not visible 
in STM) which might hinder the polymerization until bromine 
is desorbed from the surface. 

We also investigated the redistribution of the electron den¬ 
sity upon the chemisorption of the molecule on the surface. 
The difference of the total electron density between the com¬ 
bined CuTPP/Au system and the individual parts of molecule 
and Au slab is shown in Fig. 5(b),(c). The green colored 
areas correspond to an increased and the red areas to a de¬ 
creased charge density. Especially at the C atoms which form 
the C-Au bonds, the total electron density is slightly reduced, 
whereas between the Cu atom and the surface, electron den¬ 
sity increased and a small charge transfer of 0.1 electron to 
the Cu ion is found with a Bader charge analysis. This indi¬ 
cates also a direct interaction of the Cu ion in the chemisorbed 
molecule with the surface atoms due to the reduced distance 
after deformation of the molecule. 

Next, we calculated the structure of the CuTPP polymer 
in gas phase with the two coupled molecules in the periodic 
unit cell. The substrate could not be included in these cal- 
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Fig. 5 (a) Optimized geometry of debrominated CuTPPBrg on a bridge adsorption site of Au(lll). Isosurface of the difference of the total 
electron density between the combined CuTPP/Au system and the individual parts of molecule and Au slab in a side view (b) and top view (c). 
An increase of charge density is colored green, a depletion in red. 



Fig. 6 Nudged elastic band path for the polymerization reaction 
between unsaturated CuTPP molecules in gas phase and the 
dislocation of the H atom from the newly formed C-C bond to the 
unsaturated pyrrolic C atom. The change of the atomic positions at 
each step is shown and linked by arrows to the corresponding 
change of the total energy. 


culations because the polymeric network is not commensu¬ 
rable with the Au lattice and the very large simulation cell 
which would be necessary for the coincidence of both unit 
cells was not computationally feasible for us. The dimen¬ 
sions of the unit cell were optimized by comparing the total 


energy of the polymer for different cell parameters after re¬ 
laxation. The optimized size of the square unit cell was x = 
2.145nm,y = 1.485nm, which is close to the experimentally 
measured values of (2.3 ± 0.2) nm and (1.56 ± 0.20) nm. The 
calculated angle of the covalently connected pyrrole group of 
0pyr^ ~ T is reduced compared to the angle of the not cou¬ 
pled pyrrole with 0pyi- ^ 15° in gas phase, which results in the 
non-equivalent appearance for adjacent molecules in STM as 
previously discussed. 

To investigate the intramolecular displacement of a hydro¬ 
gen atom from the phenyl group to an unsaturated carbon atom 
in the pyrrole unit after the coupling, a “nudged elastic band” 
(NEB) calculation with 8 steps was performed, whereby the 
molecules in the periodic cell are already covalently coupled 
except at one phenyl and pyrrole location. The calculated re¬ 
action path is shown in Fig. 6. First in step 1 to 3, the phenyl 
group bends to get closer to the adjacent debrominated pyrrole 
group, whereby the H atom has to move out of the plane of the 
phenyl ring to enable the formation of the new covalent C-C 
bond between the unsaturated C atom and the phenyl ring. For 
this coupling reaction an energy barrier of below 0.1 eV was 
found (in gas phase). The polymerization could, thus, occur 
directly after debromination if the molecules are close enough 
and correctly oriented. Starting from NEB step 3, the H atom 
jumps on a path along the C atoms of the phenyl group to the 
still unsaturated pyrrolic C atom. The maxima of the energy 
barriers for the jumps correspond to the position of the H atom 
exactly between two adjacent C atoms. The relocation of hy¬ 
drogen to the still unsaturated C atoms leads to a large gain in 
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energy, thus, it is very likely that at least all pyrrolic C atoms 
of the CuTPP-polymer which are bent up and not in direct 
contact with the Au surface are saturated with H atoms. 


4 Conclusions 

The heat-induced coupling-limited polymerization of 
CuTPPBrg on Au(lll) by dissociation of bromine was 
studied by XPS and STM, and the results are corroborated 
with DFT calculations. For annealing until 200 °C, CuTPPBrg 
molecules debrominate and arrange in a square structure 
where the azimuthal orientation of the molecules is rotated 
by 90° in a checker board-like manner. Upon heating of 
the molecular monolayer to 300 °C or higher, C-C coupling 
between C atoms of the phenyl rings and unsaturated C in the 
pyrrole units of adjacent molecules is initiated, including the 
cleavage of a C-H bond at each phenyl group. Consequently, 
the molecules form a relatively large, covalently bonded 
network. Interestingly, shake-up peaks in XPS indicate that 
thereby the HOMO-LUMO gap opened by 0.4 eV instead of 
a decrease, which shows that the conjugation is not extended. 
The porphyrin macrocycles can, thus, still be considered as 
individual intact units although they are bonded covalently 
in a large robust 2D network. Extended two-dimensional 
polymeric sheets with high chemical and mechanical stability, 
especially from building blocks which include also magnetic 
metal ions could be promising for spintronic devices. Also, 
the lateral transport properties could be interesting for organic 
interconnections in novel molecular electronic devices. For 
this, however, it has to be possible to remove the polymer 
from the metallic surface or to decouple it, e.g. by a thin 
insulating buffer layer, which will be a big challenge. A next 
step is also to investigate the influence of the diffusion for 
the polymerization of CuTPPBrg further by choosing more 
reactive surfaces with a higher diffusion barrier, which we 
will describe in a following publication. 
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